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We investigate the instabilities that may lead to the breakdown of the Luttinger liquid in the small- 
diameter (5,0) nanotubes, paying attention to the competition between the effective interaction 
mediated by phonon- exchange and the Coulomb interaction. We use bosonization methods to 
achieve an exact treatment of the Coulomb interaction at small momentum-transfer, and apply 
next renormalization group methods to analyze the low-energy behavior of the electron system. This 
allows us to discern the growth of several response functions for charge-density-wave modulations 
and for superconducting instabilities with different order parameters. We show that, in the case of 
single nanotubes exposed to screening by external gates, the Luttinger liquid is unstable against the 
onset of a strong-coupling phase with very large charge- density- wave correlations. The temperature 
of crossover to the new phase depends crucially on the dielectric constant k of the environment, 
ranging from T c ~ 10~ 4 K (at k w 1) up to a value T c ~ 10 2 K (reached from k ~ 10). The physical 
picture is however different when we consider the case of a large array of nanotubes, in which there 
is a three-dimensional screening of the Coulomb interaction over distances much larger than the 
intertube separation. The electronic instability is then triggered by the divergence of one of the 
charge stiffnesses in the Luttinger liquid, implying a divergent compressibility and the appearance 
of a regime of phase separation into spatial regions with excess and defect of electron density. 



I. INTRODUCTION 

During the last decade, carbon nanotubes have shown 
a great potential for the development of molecular elec- 
tronics. This comes from their versatile electronic prop- 
erties, which depend on variables like the helicity of the 
tubule, the temperature, or the contacts used in the ex- 
periments. It is well-known for instance that the carbon 
nanotubes may have semiconducting or metallic behav- 
ior, depending on the roll-up direction of the wrapped 
hexagonal carbon lattice 1 . This has opened the way for 
the construction of intramolecular junctions behaving as 
electronic devices at the nanometer scaled. 

Given their reduced dimensionality, the metallic nan- 
otubes behave as strongly correlated electron systems, 
with different regimes depending on the energy scale. At 
room temperature, for instance, the strong Coulomb re- 
pulsion prevails in nanotubes of typical radius, driving 
the electron system to a state with the properties of the 
so-called Luttinger liquid, characterized by the absence 
of electron quasiparticles at the Fermi level. This feature 
manifests in the power-law dependence of observables 
such as the tunnelling density of states, whose suppres- 
sion at low energies has been actually observed in mea- 
surements of the conductance in the carbon nanotubes 2 ^. 

At low temperatures, the behavior of the carbon nan- 
otubes depends on the contacts used in the transport 
measurements. When the contacts are not highly trans- 
parent, there is in general a suppression of the zero-bias 
conductance and the differential conductivitiy, charac- 



teristic of the Coulomb blockade regime. On the other 
hand, there have been experiments using high-quality 
contacts where it has been possible to observe super- 
conducting correlations in the carbon nanotubes^. The 
most remarkable signature has been the appearance of 
supercurrents in nanotube ropes suspended between su- 
perconducting electrodes^. This feature may be under- 
stood as a consequence of the proximity effect, by which 
Cooper pairs are formed in the nanotubes near the su- 
perconducting electrodes^. Moreover, superconducting 
transitions have been measured in ropes suspended be- 
tween metallic, nonsuperconducting electrodes^. These 
observations imply the existence of a regime with a rele- 
vant attractive component in the electron-electron inter- 
action, coming presumably from the exchange of phonons 
between electronic currents^. 

^From a theoretical point of view, it is well-known 
that the Luttinger liquid regime may break down in the 
carbon nanotubes, due to a variety of low-temperature 
instabilitie»i2*ii*i2ii2iii When the effective phonon- 
mediated interactions are taken into account, it can be 
shown that the superconducting correlations may grow 
large upon suitable screening of the Coulomb interac- 
tion in the nanotubes^. This has raised the hopes that 
such correlations could be amplified by considering nan- 
otube geometries with enhanced electron-phonon cou- 
plings. The case of the small-diameter nanotubes, with 
radius down to 0.2 nm, is particularly interesting in that 
respect, since the large curvature of the tubule increases 
significantly the coupling to the elastic modes of the nan- 
otube latticeiSiii. 

The transport properties of small-diameter nanotubes 
have been studied in the experiment reported in Ref. ITsl 
where it has been claimed evidence for a superconduct- 
ing transition at about 15 K. It is uncertain however 
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FIG. 1: Representation of the Brillouin zone of the (5,0) nan- 
otubes showing the position of the Fermi points belonging to 
the low-energy degenerate subbands (denoted by 1,2) and to 
the zero angular-momentum subband (denoted by 0). 



whether the observed features should be ascribed to the 
presence of nanotubes with the (5,0) or the (3,3) geom- 
etry, both being consistent with the measured nanotube 
radius. Moreover, it remains also unclear the physical 
meaning of the estimated transition temperature, as the 
nanotubes in the experiment are placed in the channels 
of a zeolite matrix that does not allow coherent electron 
hopping between the nanotubes. 

The aim of the present paper is to investigate the insta- 
bilites that may lead to the breakdown of the Luttinger 
liquid in the small-diameter zig-zag nanotubes. We are 
going to focus on the analysis of the (5,0) nanotubes, 
which deserve special attention as they have three sub- 
bands crossing the Fermi level (in the undoped system). 
The position of the different Fermi points in momentum 
space is shown in Fig. ^ They correspond to two de- 
generate subbands with opposite values of the angular 
momentum around the nanotube axis (denoted by 1,2) 
and a third subband with zero angular momentum (la- 
belled by in the figure). The Fermi velocity is different 
in the nondegenerate and the degenerate subbands, be- 
ing smaller in either case than for nanotubes of normal 
radius and leading to an enhanced density of states at 
low energies. For the degenerate subbands, the Fermi 
velocity becomes vf ~ 2.8 x 10 5 ms" 1 , while for the non- 
degenerate subband we have z4 0) « 6.9 x 10 5 ms . A 
diagram with the different shapes of the subbands along 
the momentum in the longitudinal direction is shown in 
Fig. El 

Together with the enhanced density of states at low en- 
ergies, the large curvature of the tubule makes the (5,0) 
geometry the most appealing instance to study the effects 
of large electronic correlations in the carbon nanotubes. 
We will apply renormalization group methods^ to dis- 
cern the low-energy instabilities of the electron system 
with three subbands at the Fermi level, extending the 
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FIG. 2: Schematic representation of the dispersion of the low- 
energy degenerate subbands (superposed in the curve labelled 
by 1 and 2) and the zero angular-momentum subband (la- 
belled by 0) of the (5,0) nanotubes. The Fermi level corre- 
sponds to e — (in the undoped system). 

analysis carried out in Refs. and|2(] for typical nan- 
otubes with two low-energy subbands. This will allow us 
to study the growth of several response functions marking 
the tendency towards long-range order, including charge- 
density-wave (CDW) modulations with different values 
of the momentum and superconducting instabilities with 
different order parameters. 

Our analysis will stress the relevance that the screen- 
ing of the Coulomb interaction has for the development 
of different instabilities in the small-diameter nanotubes. 
For a single nanotube, placed in general above some sub- 
strate, there is little reduction of the Coulomb potential, 
which remains long-ranged in the one-dimensional (ID) 
systemSi. The Coulomb interaction has then a strong 
effective coupling, much larger than that of the effec- 
tive phonon-mediated interaction, so that the Coulomb 
repulsion plays the dominant role dictating the elec- 
tronic properties about room temperatures. We will use 
bosonization methods^ to achieve an exact treatment of 
the contributions from the Coulomb interaction at small 
momentum-transfer, capturing the effects of the strong 
Coulomb repulsion at the nonperturbative level. The 
contributions from the effective phonon-exchange inter- 
actions will be analyzed next, assessing their enhance- 
ment at low energies from the nontrivial scaling of the 
backscattering interaction channels. 

On the other hand, we will pay also attention to the 
particular conditions in the experiment reported in Ref. 
Ha . Each nanotube is embedded there in a large three- 
dimensional (3D) array of nanotubes. Thus, the elec- 
trostatic coupling between charges in the different nan- 
otubes leads to a large screening of the Coulomb poten- 
tial. We will study this effect by using a kind of general- 
ized RPA scheme, taking into account the Coulomb inter- 
action between all the nanotubes in the 3D array. We will 
show that, over distances much larger than the intertube 
separation, there is a regime where the nanotube array 
screens effectively as a 3D system, rendering short-ranged 
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the Coulomb potential within each nanotube. This is the 
most important effect of the environment in the partic- 
ular experimental setup described in Ref. 0. We will 
see that the strong reduction of the Coulomb interaction 
leads then to the prevalence of charge instabilities with 
similar character to the Wentzel-Bardeen singularity^. 

We may compare the results of our investigation with 
those obtained for the small-diameter zigzag nanotubes 
by means of ab initio simulations 2 ^ and mean-field-like 
calculations^. In Ref. |2^ a superconducting instability 
has been found to be dominant in the (5,0) nanotubes, at 
a scale of about 1 K. On the other hand, a CDW insta- 
bility has been identified in the same geometry at room 
temperature in Ref. 0- Our results stress that the scale 
and the type of electronic instability are very sensitive to 
the degree of screening applied to the nanotubes. We will 
find indeed that single (5,0) nanotubes (with just some 
screening from nearby gates) undergo the breakdown of 
the Luttinger liquid regime in favor of a CDW modu- 
lation at low energies. This is quite different from the 
instability found in the 3D array of nanotubes. In either 
case, we will see that the onset of the electronic instability 
takes place at temperatures ranging from T c ~ 10~ 4 K 
up to ~ 10 2 K, depending on the dielectric constant of 
the environment. 



II. LUTTINGER LIQUID APPROACH TO 
ELECTRONIC PROPERTIES 

As a first approximation to the electronic properties of 
the small-diameter nanotubes, we begin our analysis by 
focusing on the effects of the strong Coulomb repulsion. 
This is the dominant interaction in scatteri ng p rocesses 
at low momentum-transfer. Following Ref. [TJ, we take 
the Coulomb potential in the wrapped geometry as 

Vc(r r') = e2/K 

yj(x- x') 2 + 4i? 2 sin 2 [(y - y') 2 /2R 2 } + a 2 

where a z ~ 1.6 Aand R is the nanotube radius. The effec- 
tive strength of the interaction depends on the dielectric 
properties of the environment hosting the nanotubes. In 
general, we will encode the screening effects from exter- 
nal gates in the dielectric constant k. We have to bare 
in mind, however, that this procedure is suitable for the 
description of single nanotubes, while it has to be ap- 
propriately improved for systems made of a manifold of 
nanotubes, as we will see in Sec. IV. 

The Fourier transform Vc(k, q) of the potential (JJJ re- 
flects the long-range character of the interaction at small 
momentum-transfer. Thus, it has a logarithmic depen- 
dence at small longitudinal momentum k, which is char- 
acteristic of the ID Coulomb potential in momentum 
spaced: 

Vb(fc,0)^log(^), (2) 



k c is in general of the order of the inverse of the nan- 
otube radius R, as it is the memory that the electron 
system keeps of the finite transverse size, after projection 
of the 3D potential onto the longitudinal direction of the 
tubule. The strength of the Coulomb interaction can be 
estimated from the dimensionless ratio e 2 / kvf, which for 
the small-diameter nanotubes is of the order of ~ 8/k. 
This gives a measure of the relevance of the Coulomb in- 
teraction, which turns out to be in the strong-coupling 
regime in the forward-scattering channels. 

Before going ahead, we proceed next to make a com- 
plete catalogue of the different interaction channels. We 
first borrow the classification that has been already made 
for metallic nanotubes with a pair of subbands at the 

Fermi \eve&. Thus, we introduce respective coupling 

(7) 

constants g\ for the channels involving processes be- 
tween the two inner subbands of the small-diameter nan- 
otubes. The lower index discerns whether the interacting 
particles shift from one subband to the other (i = 1), 
remain at different subbands (i = 2), or they interact 
within the same subband (i — 4). The upper label fol- 
lows the same rule to classify the different combinations 
of left-movers and right-movers. 

The above set of couplings gf has to be supplemented 
with the couplings for the channels involving processes 
between a particle at one of the degenerate subbands and 
another at the nondegenerate subband. Now we can dis- 
cern between interactions that keep each particle within 
its respective subband, to which we assign a coupling 
/( 2 )(/( 4 )) for particles of opposite (like) chirality, and 
backscattering interactions that lead to the exchange of 
the subbands for the two particles, which we label with 
a coupling /W. Moreover, we have also the channels in 
which the particles interact within the outer subband, 
and that we will label as in the usual <?-ology descrip- 
tion, in terms of the couplings g^\g^ and g^ . Finally, 
there are also interaction processes in which two particles 
with about zero total momentum exchange their position 
from the outer to the inner subbands, or vice versa, as 
depicted in Fig. |3J We will assign the coupling uf to 
that kind of interaction with no change of chirality of the 
particles, and the coupling ug when there is a change 
in their chirality. It is worthwhile to remark that, de- 
spite the role they may have in the enhancement of some 
type of correlations, these interactions labelled by uf and 
ub have been overlooked in previous studies of the (5,0) 
nanotubes 2 ^. 

Focusing again on the forward-scattering interactions, 
we may represent their contribution to the hamiltonian 
of the electron system introducing the charge and spin 
density operators 

Pn(x) = ^(9l ir (x)9 ri ^x)+9l ii (x)9 riL {x)) (3) 

<r*(x) = -^(*t n (x)* ri1 (x) - ^(x^ri^x)) (4) 

which correspond to the different electron fields ^ r ia for 
the linear branches about the Fermi points shown in Fig. 
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FIG. 3: Diagrams showing the shift in momentum space of 
the interacting electrons and the momentum-transfer involved 
in different backscattering processes in the (5,0) nanotubes. 



13 We adopt a notation in which the index r = L,R is 
used to label the left- or right-moving character of the 
linear branch, and the index i = 0, 1, 2 to label the sub- 
band. We will assume that the interaction mediated by 
phonon exchange (as well as the Coulomb interaction) 
does not depend on the spin of the interacting electrons, 
so that we will carry out the discussion in terms of the 
charge density operators. 

For the sake of organizing the different forward- 
scattering interactions, it is convenient to define the sym- 
metric and antisymmetric combinations of corresponding 
density operators in the two degenerate subbands: 

PR±(k) = -j={p R i{k) ± pm{k)) (5) 

PL±{k) = ^=( PL2 (k) ± PL1 (k)) . (6) 

With this change of variables, the hamiltonian for the 
forward-scattering interactions can be written in the form 

^ /"fee 

HfS = rt'F / dk 2J ^ Pri(k)pri(-k) 



,(0) 



-fee 



r=L,Ri=± 

dk 2J Pro(k) Pro(-k) 



r=L,R 



+ \ f C k ^2( PR+ (k) (g^+g^)p R+ (-k) 
+p L+ (k) ( 3 f + PL+ (-k) 



+PR-W (sf ~ 92 ] ) PR-(-k) 
+p L .(k) (g^ - g ( 2 i} ) PL-(-k) 
+2 PR+ (k) (g^+g^PL+i-k) 
+2p R .(k) (g™ - 5 f) p L -(-k)) 

+ \ ^2(p R0 (k)g^ pm{ -k) 
2 J-k c 27r v 

+ PLa (k) g^p L0 (-k) + 2p m (k) gWp L0 (-kj) 

+ \ f C k ^*V2(pR + (k) /W Pm (-k) 

+p L+ (k) f ( ^p L0 (-k)+p R+ (k) fWp L o(-k) 
+ PL+ {k) fW PR0 (-k)) , (7) 

where k c stands again for the momentum cutoff dictated 
by the transverse size of the nanotube. 

In order to diagonalize the hamiltonian lj?Jl. we intro- 
duce the fields $+(x), $_(x), $o(x) (and their respec- 
tive conjugate momenta, 11+ (a;), II„ (x), IIo(x)) having 
the following relation with the electron density opera- 
tors: 

d x $+(x) = y/^(p L+ (x) + p R+ (x)) (8) 
d x $>-{x) = V^(pl-(x) + pr-{x)) (9) 
d x ® {x) = Vn(p L o(x) +prq(x)). (10) 

In terms of these fields, we can write the hamiltonian Q 
in the form 

Hps = \u+ f dx{K + (U + (x)) 2 + (l/K + )(d x <S> + (x)) 2 ) 



dx (A_(LF_(x)) 2 + (l/^-)(^$_(x)) 2 ) 



+ -u I dx (A (n (x)) 2 + {l/K ){d x ^ {x)f) 

i r , 2V2 



+ ^j dx -^r ( n+(x) (/ (4) -/ (2) ) n oW 

+d x $+(x) (/ (4) +/ (2) )^$ (x) 



(11) 



The renormalized velocities u + ,u-,uq and charge stiff- 
nesses K + , K- , Kq are given by the equations 



u±K± 


— VF - 


u±/K ± 


— Vp - 


UqKq 


(0) 

= vy 




(0) 

= v F 



(14) 
(15) 



At this stage, the hamiltonian l|llf> can be brought to 
diagonal form by: i) applying a simple canonical trans- 
formation 



= \/A + <i> i 



n. 



1 
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$o = v^o^o , n 



=n 



(16) 



and ii) making an additional rotation to disentangle the 
(+,0) sector: 



$o = 

n = 

$+ = 
n+ = 



n 



•Sf) 







IL 



VP 



n 



+ 



n„ 



(17) 



where the parameters cq, Sq, c + , s+, a,/3, /i, v are evalu- 
ated in Appendix A. The important point is that, at 
the end, the excitations of the system governed by the 
hamiltonian (|ll(l are given by charge fluctuations, with 
velocities strongly renormalized by the interactions, and 
which constitute the whole spectrum together with the 
spin fluctuations propagating with the unrenormalized 
velocities Vp and Vp . 

The picture that we have just developed corresponds 
to the Luttinger liquid regime of the electron system, 
in which the main physical properties are given by the 
charge stiffnesses in the different sectors and by the veloc- 
ities of the charge and spin fluctuations^. These param- 
eters may be used, for instance, to characterize the com- 
pressibilities in the symmetric, antisymmetric and zero 
angular-momentum sectors, which turn out to be pro- 
portional to the respective charge stiffnesses. 

The Luttinger liquid behavior in the carbon nanotubes 
holds at energies where there is approximate linear dis- 
persion of the different subbands crossing the Fermi level. 
This places typically an upper scale E c of the order of 
~ 0.1 eV. In the case of single nanotubes, the domi- 
nant contribution to the forward-scattering couplings in 
(Q) comes from the strong Coulomb potential in © , with 
some infrared cutoff which is dictated in general by the 
longitudinal size of the nanotubeii. This leads to values 
of the charge stiffnesses that are typically smaller than 1 
and in agreement with the experimental observations of 
Luttinger liquid behavior in the carbon nanotubes 2 *^. For 
energy scales much lower than E c , however, the rest of 
interaction channels neglected in the above analysis may 
become as important as the forward-scattering channels. 
This is possible as the strength of the ID interactions 
depends in general on the energy scale of the interaction 
processes^. Electronic instabilities are then expected 
when some of the backscattering interactions leave the 
weak-coupling regime as the temperature or other rele- 
vant energy scale is lowered, as we analyze next in the 
case of the small-diameter (5,0) nanotubes. 



III. BACKSCATTERING INTERACTIONS AND 
LOW-ENERGY INSTABILITIES 

The small-diameter nanotubes are systems where the 
backscattering interactions may lead to a significant en- 
hancement of the electronic instabilities. Due to the 
large curvature of the tubes, the electron-phonon cou- 
plings turn out to be much larger than for nanotubes 
of typical radius. This leads consequently to a larger 
electron-electron interaction mediated by the exchange 
of phonons. In the case of single nanotubes, this ef- 
fective interaction is anyhow overcome by the strong 
Coulomb repulsion given by the potential in the 
forward-scattering channels. The opposite situation is 
found however in the backscattering channels, as we dis- 
cuss in what follows. 

In general, the exchange of phonons between electronic 
currents gives rise to a retarded interaction, which can 
be represented by the effective potential obtained by in- 
tegrating out the phonons in the many-body theory^. If 
we denote the electron-phonon couplings by gw(k) (la- 
belling by the subbands of the incoming and the out- 
going electron), we can express the potential for the ef- 
fective electron-electron interaction in the form 



V; 



(k,u) = -2ffii/(fc)5jj-/(-fe)- 



(18) 



where LUk stands for the energy of the exchanged phonon 
with momentum k. The typical energy of the phonons in 
the optical branches (or in the acoustic branches at large 
momentum-transfer) is of the order of ~ 0.1 eV 25 , and 
therefore comparable to the energy cutoff E c of the ID 
electron system. This allows us to take the interaction 
arising from Ijl8|l as a source of attraction, in the energy 
range where the ID model makes sense. 

The electron-phonon couplings in the (5,0) nanotubes 
have been analyzed with great detail in Ref. |2^. It 
has been shown there that the couplings with greater 
strength correspond to intraband processes, for which 
the momentum-transfer is in the longitudinal direc- 
tion. Following for instance the results of that refer- 
ence for the backscattering processes within subband 1 
(or subband 2), we find that the strength of the effec- 
tive electron-electron interaction with 2k p momentum- 
transfer is given by the coupling 



(gtt&kp)) 2 /v F 4 



0.6 



(19) 



after summing over the different phonon modes with 
nonvanishing electron-phonon coupling 25 . Similarly, for 
backscattering processes within the nondegenerate sub- 
band we find the effective coupling 



^W)) / 



W) 

VFLU 2k^ 



1.1 



(20) 



The first of these values represents the contribution of 
the attractive phonon-mediated interaction to the cou- 
plings an d <?4 j while the second value provides the 
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contribution to the coupling Making similar esti- 

mates of the electron-phonon couplings and frequencies 
at large transverse momentum-transfer, we have found 
that the couplings for the rest of backscattering interac- 
tions g± ,gx , f^\uF and ub are given by w Q.'bvp. 

To each of these backscattering contributions, one has 
to subtract an opposite correction from the Coulomb 
interaction. This is given by the Fourier transform 
Vc(k,q) of the potential in at the corresponding 
momentum-transfer, with an additional factor of reduc- 
tion that comes from the particular structure of the Bloch 
functions in the small-diameter nanotubes 25 . In gen- 
eral, the larger contributions from the Coulomb potential 
are found in the channels with the smaller momentum- 
transfer. The forward-scattering channels deserve special 
consideration, as the potential © has to be evaluated 
with a cutoff at small k. As already mentioned, this is 
provided in general by the length L of the experimen- 
tal samples, so that a sensible estimate is given typi- 
cally by a spatial average of the potential Vc(k «0,0)a 
(2e 2 /K)log{k c /k ), with k ~ 1/L ~ l(T 3 fc c . 

In this section we analyze the instabilities of single 
nanotubes, so that we can disregard for the moment 
any kind of intertube interaction. Then, taking into ac- 
count the contributions from the phonon-exchange and 
the Coulomb interaction, we obtain the following values 
for the couplings in the different interaction channels: 

aflvp = g^/vp = g^/vp = g^/vF « 100/* (21) 

g^/vp = g {2) /v F = f {4) /v F = f (2) /v F « 100/k (22) 

g [ l ] /v F = g^/vp ~ -0.6 + 3.6/ k (23) 

gW/v F » -1.1 + 1.1/k (24) 

f^/vp = u B /v F « -0.5 + 0.5/k (25) 

up/vp ~ -0.5 + 0.3/k (26) 

g^/vp = g ( i ] /vF ~ -0.5 + 0.05/k .(27) 

We note again that the different contributions from the 
Coulomb interaction (with a 1/k reduction in each case) 
correspond to the pertinent choices of the momentum- 
transfer in the evaluation of Vc(k, q). 

The relevance of the backscattering interactions comes 
from the fact that they give rise to quantum corrections 
that depend logarithmically on the energy scale. In the 
case of the g^ interaction, for instance, the bare cou- 
pling is corrected to second order by the diagrams shown 
in Fig. 21 By following the renormalization group pro- 
gram, we translate the logarithmic dependence on en- 
ergy of the terms in the diagrammatic expansion into the 
scale-dependence of renormalized coupling constants 19 . 
Up to terms quadratic in the backscattering interactions, 
we find the complete set of scaling equations 



+{(3/2){u 2 F + u 2 B )) 



(29) 



dg\ 



(i) 
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a (2 
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dl 



(Si Si +Si S2 i pu F u 
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dgW 
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dl 

duF 

~dT 
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di 



(1-^)S^-— (2^-si 1 ^ 



+Si 2) Si 1) + Pufub) 



(30) 



~^ 1) +^gjP+M (3D 



TTVF 

-^—(a (1) a {1) ~a {2) a {2) ) 
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TTVF 

- — {(ig [1) g (1) + 2u f ub) 
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(/ (1) / (1) -4) 



2ttvf 

a 
irvp 

Auf - 



(f {1) f {1) +U 2 B-UFU B ) 

' I (2) 
(Si UF 



2ttvf 



+ (S^ ) +S[ 1) +/3S (1) W) 



/ (2) 
(Si UB 



Au B - - 

Zttv p 

+(gi 1] + g{ 1} + Pg (1) >F) 

a -f^{u F -2u B ), 



TTVF 



(32) 
(33) 
(34) 
(35) 
(36) 
(37) 

(38) 
(39) 



(1 



\ (2) 

)si 



TTVF 



where (3 = vf/v f , a — 2/(1 + v F /vf) and the anoma- 
lous dimension A depends on K + , K-, Kq and as 
shown in Appendix A. In the above equations, the vari- 
able I stands for minus the logarithm of the energy (tem- 
perature) scale measured in units of the high-energy scale 
E c of the ID model (of the order of ~ 0.1 eV). 

The equations for the g^ couplings correspond to 
those already written for a model with two subbands in 
Ref. l20i Here we have incorporated a nonperturbative 
improvement of the equations by writing the exact de- 
pendence of the anomalous dimensions on the forward- 
scattering couplings through the K + , K- and Kq param- 
eters. This is essential to deal with the strong-coupling 
regime of the forward-scattering interactions in the case 
of single small-diameter nanotubes. Another important 
novelty is that the equations for the couplings in the de- 
generate and the nondegenerate subband mix through 
the couplings uf and ub, which had been overlooked 
however in previous analyses of the (5,0) nanotubes. 

In order to determine the electronic instabilities that 
may appear in single small-diameter nanotubes, we have 
solved the set of scaling equations (|28|I - H39|I . taking ini- 
tial values for the couplings according to the above dis- 
cussion. The essential point is that the scaling equations 
reflect the strong screening that the couplings undergo 
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FIG. 5: Plot of xcdw{}c) as a function of re. Here l c = 
\og(E c /u c ) is the scale at which the flow breaks down. The 
three regions with the corresponding divergent Luttinger liq- 
uid parameter are also indicated. 



FIG. 4: Second-order diagrams with logarithmic dependence 
on the frequency renormalizing the interaction. The full 
(dashed) lines represent the propagation of electrons with 
right (left) chirality, and the wavy lines stand for the inter- 
actions. The labels 0,1,2 denote the respective low-energy 
subbands of the (5,0) nanotubes. 

at low energies, which leads them to a regime of large at- 
traction. To find the dominant instability in the electron 
system, one has to look at the behavior of the different 
response functions x, which are catalogued in Appendix 
B. The regime of attraction leads to the large growth of 
some of the x> this points at a tendency towards long- 
range order, that cannot be completed anyhow at any 
finite frequency in the ID system. The scaling of the 
interactions is cut off at some low-energy scale co c for 
which one of the charge stiffnesses j3, fi or K- diverges. 
We characterize the dominant instability by identifying 
the response function that reaches the largest value at 
the scale u> c . 

By varying the dielectric constant, we have found two 
different regimes, namely a regime characterized by the 
divergence of (for 1 < n < 1.6 and k > 3.4), and a 
regime with the divergence j3 (for 1.6 < k < 3.4). The 
largest growth for the respone functions corresponds to 
Xcdw i with large momentum connecting opposite Fermi 
points of the degenerate subbands (see Appendix B). We 
can understand this finding by looking at the compe- 
tition between the phonon-exchange and Coulomb con- 
tributions in the different channels. At small k, the cou- 
plings for large transverse momentum-transfer (g^\gi) 
are attractive from the very beginning (see Ea. (|27|l 'l. In 

particular, the flow is dominated by g± , which is driven 
to large attraction, leading to a rapid growth of xcdw- 
Anyway, as shown in Fig|51 the xcdw grow large only 
for small values of k (k ^ 1.4-^2 for nanotubes in typical 
conditions). At large k the maximum value for the all 
the respone functions remain close to 10, signaling that 




2L , , , 1 

2.5 5 7.5 10 12.5 15 17.5 20 
K 

FIG. 6: Plot of minus the logarithm of the energy scale l c — 
log(E c /Lu c ) at which the flow breaks down, as a function of k. 



the tendency to ordering is weak. 

In FigEl we have plotted the critical value l c = 
\og(E c /u! c ) at which the flow breaks down as a function 
of k. It appears that the energy scale lo c becomes quite 
sensitive to the value of the dielectric constant, implying 
that the onset of the electronic instability can be found 
at temperatures ranging from T c ~ 10~ 4 K (at k w 1) 
up to a value T c ~ 10 2 K (reached from n w 10). Fi- 
nally, we observe from the comparison in Fig. E| that, 
while some charge and spin-density-wave response func- 
tions grow large by approaching the critical value l c (at 
least for small k), the superconducting response func- 
tions for different order parameters remain all small. This 
finding seems to rule out the possibility of having super- 
conducting correlations in the (5,0) nanotubes, at least 
under the physical conditions considered in the present 
section. 
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IV. INTERTUBE SCREENING AND 
INSTABILITIES IN A 3D ARRAY OF 
NANOTUBES 



While the above framework provides an accurate de- 
scription of the instabilities of single small-diameter nan- 
otubes, the analysis has to be conveniently modified when 
the electron system is composed of a manifold of nan- 
otubes. This is the situation in the experiments reported 
in Ref. where the samples contain large 3D arrays 
of small-diameter nanotubes. Systems that are formed 
by the assembly of a large amount of nanotubes, as is 
also the case of nanotube ropes, require in general an ad- 
ditional analysis of the interactions arising between dif- 
ferent nanotubes. We will now focus on the conditions 
of the experimental samples studied in Ref. ITsl which 
present a negligible intertube tunneling but have a sig- 
nificant coupling between electronic currents in different 
nanotubes. 

Specifically, we will consider the case of an array of 
small-diameter nanotubes, arranged as a triangular lat- 
tice as viewed from a cross-section of the 3D array. We 
will denote by d the distance between nearest-neighbor 
nanotubes in such a lattice, having in mind that, for 
the samples studied in Ref. fl8l its value is given by 
d w 1 run. The coupling between different nanotubes 
at positions 1 and 1' (measured in the cross-section of 
the array) comes from the Coulomb interaction, due to 
its long-range character. It may couple currents with 
large transverse separation |1 — l'| as long as the longitu- 
dinal momentum-transfer k becomes as small as |1— 1'| . 
The Coulomb potential Vij' (fc) between currents in differ- 
ent nanotubes 1 and 1' is actually obtained by partially 
Fourier-transforming the 3D Coulomb potential in the 



FIG. 8: Second-order process with logarithmic dependence on 
the frequency renormalizing intraband interactions at a given 
nanotube 1 through the coupling with the nearest-neighbors 
1' in a 3D array of nanotubes. The dashed lines (without 
arrow) stand for interactions between electronic currents in 
nearest-neighbor nanotubes. 



longitudinal direction: 



2e 2 



V h y(k)^—K (\l-l'\k), 
K 



(40) 



Kq(x) being the modified Bessel function, which diverges 
logarithmically as x — > and is exponentially suppressed 
as for x > 1. 

It is clear then that the nanotubes in the array 
may screen efficiently the forward-scattering interac- 
tions within each nanotube. This effect can be studied 
by extending the Random Phase Approximation (RPA) 
scheme to incorporate the electrostatic coupling between 
all the nanotubes in the array, as described in Ap- 
pendix C. The main conclusion of this study is that 
the screened Coulomb potential V l j (k) within each nan- 
otube becomes finite in the limit of vanishing momen- 
tum k —>■ 0, reaching a saturation value VJj (k = 0) « 
0.08e 2 . This corresponds to a dimcnsionless coupling 
Vjj (k — 0)/vp 0.65, which is about two orders of mag- 
nitude smaller than the one for single nanotubes. In this 
case the phonon-exchange contribution competes with 
the Coulomb repulsion also in the intratube forward- 
scattering channels, changing qualitatively the physical 
picture. 

One more important effect is that the coupling be- 
tween currents in neighboring nanotubes may also cor- 
rect the intratube backscattering interactions with small 
momentum-transfer. This is due to the fact that, for a 
longitudinal momentum-transfer 2kp, the Coulomb po- 
tential between nearest-neighbor nanotubes (such that 
|1 — 1'| = d) still has a nonnegligible strength, V\j'(2kp) ~ 
0.037e 2 . The important point is that this interaction 
gives rise to processes of the type depicted in Fig. |S1 
which depend on the energy scale and therefore intro- 
duce important renormalization effects for the intratube 
interactions at low energies. As observed from Fig. the 
scaling analysis requires the definition of new backscat- 
tering interactions that couple electron currents in differ- 
ent nanotubes. We will assign the couplings and 
to the new channels represented in Figs. EI a )-(b). As 
we are going to see next, these new couplings mix upon 
renormalization with two more channels corresponding 
to forward-scattering interactions between currents with 
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FIG. 9: Intertube interactions arising from the coupling be- 
tween electronic currents in nearest-neighbor nanotubes 1 and 
1' of a 3D array. The meaning of the different lines is the same 
as m Figs. Handel 



different chirality in nearest-neighbor nanotubes. These 

(2) (2) 

will be labelled by the couplings g\ and g% , as shown 
in Figs. Etc)-(d). 

In order to determine the perturbative scaling of the 
new intertube interactions g^ and g£ , we pay attention 
to the second-order corrections that have logarithmic de- 
pendence on the energy scale. These have to be made of 
particle-particle processes or particle-hole loops with the 
nesting momenta 2Uf or 2k p . The different diagrams 
are shown in Fig. ^| We see that some of them involve 

~(2) ~(2) 

the new intertube interactions g\ and g\ . A consis- 
tent renormalization demands then the analysis of the 
scaling of this pair of couplings. This is dictated again 
by diagrams that depend logarithmically on the energy 
scale, depicted in Fig. ^Jto second order in perturbation 
theory. 

The main role of the intertube backscattering interac- 
tions g± and g^ is to screen the intratube backscatter- 
ing interactions through processes of the type shown in 
Fig. [S] Thus, the scaling equations for g^ 1 and gjp get 
new contributions in the 3D array of nanotubes, taking 
now the form 



dl 



dl 



1 



K- 7TV F 



1 



(1) (2) 

94 9i 



+g{ 2) g[ 1) +Pu F u B + l2gi 1) g { 2 1) ) (41) 



irvp 



(2) (1) 

9i 92 



(42) 



We observe that the new terms are enhanced by a fac- 
tor proportional to the number of nearest-neighbors of 
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FIG. 10: Second-order diagrams with logarithmic dependence 
on the frequency renormalizing the intertube interac- 
tion. The wavy lines stand for intratube interactions and the 
dashed lines (without arrow) for interactions between nearest- 
neighbor nanotubes 1, 1' in a 3D array. 



each nanotube in the 3D array. The scaling equations 
for the new intertube couplings follow from the diagrams 
depicted in Figs. ITUlandlTTl 
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(46) 



Let us now discuss the initial conditions for the new 
set of scaling equations. The intratube backscattering 
couplings have the same initial values as in the preceding 
section since, according to the above discussion, the inter- 
tube screening in these channels is already incorporated 
in the scaling equations. Conversely, in the forward- 
scattering channels, the screening may be represented by 
finite diagrammatic corrections, which can be summed 
up with the RPA approach described in Appendix C. As 
mentioned before, the Coulomb contribution is given by 
the dimensionless coupling V^[\k = 0)/vp ~ 0.65. The 
phonon-exchange contribution at vanishing momentum- 
transfer can be estimated to be very similar to the value 
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FIG. 11: Second-order diagrams with logarithmic dependence 
on the frequency renormalizing the intertube g^ 1 interaction. 
The dashed lines (without arrow) stand for interactions be- 
tween nearest-neighbor nanotubes 1, 1' in a 3D array. 



found from (|19fl , showing that there is now a substantial 
suppression of the effects of the Coulomb repulsion in the 
intratube forward-scattering channels. 

On the other hand, the intertube forward-scattering 
interactions are also screened due to the electrostatic cou- 
pling among the nanotubes in the array. These screening 
effects can be evaluated again within the RPA scheme 
in Appendix C, with the result that they render finite 
the intertube screened Coulomb potential at k — > 0. 
Furthermore, the values of the intertube backscattering 
couplings g^ and gp can be obtained from l|40|) with 
|1 — l'| = d. The new initial conditions obtained in this 
way read: 



gf ] /v F = gf ] /v F ~ 0.002k 
gW/vir = g<P/v F « 0.28/re 



(47) 
(48) 



We observe that the re-dependence of the screened 
forward-scattering couplings is qualitatively different 
with respect to the nonscreened interactions. In particu- 
lar, we find that the intertube couplings have an approx- 
imate linear dependence on re, at least up to re ^ 5 (see 
Appendix C). 

The numerical integration of the scaling equations 
shows now a new physical regime, where the instability 
of the system is not characterized by the large growth of 
any of the response functions, but only by the divergence 
of the charge stiffness j3. We have plotted in Fig. ^] a 
typical flow of the dominant response functions up to the 
energy scale at which the divergence of the charge stiff- 
ness occurs. We note that the screening effects produced 




FIG. 13: Plot of minus the logarithm of the energy scale l c = 
\og(E c /uj c ) at which the flow breaks down, as a function of re. 
The divergent Luttinger liquid parameter is also indicated. 



by the environment of nanotubes suppress the tendency 
to CDW ordering, and enhance the superconducting cor- 
relations. Anyway, the values of Xs's ^° n0 ^ S row large , 
making very unlikely that they may give rise to any ob- 
servable feature. 

Otherwise, the divergence of (3 signals the onset of a 
regime of strong attraction, as it implies the divergence 
of the compressibility and density-density correlator in 
the corresponding sector and the vanishing of the renor- 
malized Fermi velocity u + / (3 as well. This points at 
the development of phase separation into spatial regions 
with different electronic density, in close analogy with the 
physical interpretation of the Wentzel-Bardeen singular- 
ity. The critical energy scale for this instability is shown 
in Fig. For the zeolite matrix of Ref. 0, the estimate 
of the dielectric constant gives re w 2 ~ 4, which corre- 
sponds to a transition temperature T c sa 3 K 4- 20 K, in 
qualitative agreement with the value T c sa 15 K observed 
experimentally. 
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V. CONCLUSION 

In this paper we have studied the electronic instabil- 
ities of the small-diameter (5,0) nanotubes, by analyz- 
ing the competition between the Coulomb interaction 
and the effective interaction mediated by the exchange 
of phonons. We have built our framework on the ba- 
sis that the Luttingcr liquid is what characterizes the 
normal state of the carbon nanotubes. This is a key as- 
sumption, since the low-energy excitations of the ID elec- 
tron system are given in the normal state by charge and 
spin fluctuations, with absence of electron quasiparticles 
in the spectrum. This is consistent with the power-law 
suppression of the conductance and the differential con- 
ductivity that has been measured in several transport 
experiments on carbon nanotubes^*^. The deviations 
from this picture have to be understood as perturbations 
of the Luttinger liquid behavior, which in general are 
nominally small as the strong Coulomb potential at small 
momentum-transfer provides the dominant interaction in 
the system. 

There are however instances in which the enhance- 
ment of the backscattering interactions (corresponding 
to some nesting momentum between Fermi points) leads 
to the breakdown of the Luttinger liquid picture. This 
certainly happens at sufficiently small temperature, since 
the backscattering interactions are amplified at low ener- 
gies until they reach the strong-coupling regime 19 . There 
are also experimental conditions in which the strong 
Coulomb repulsion is largely screened. The screening 
effects are particularly important when the nanotubes 
form large assemblies 2 ^. This should explain for instance 
the appearance of a regime with superconducting corre- 
lations at low temperature in massive ropes. In this con- 
text, the increase in the curvature of the small-diameter 
nanotubes has to give rise to larger electron-phonon 
couplings 16 , implying in turn a natural enhancement 
of the backscattering interactions mediated by phonon- 
exchange. Thus, there have been good prospects to 
observe the effects of large electron correlations in the 
small-diameter nanotubes, specially in the (5,0) geome- 
try which leads to a relatively large density of states at 
the Fermi level. 

We have shown that the level of screening of the 
Coulomb interaction plays an important role in the de- 
velopment of the low-energy instabilities in the (5,0) nan- 
otubes. In the case of single nanotubes exposed to the 
screening by external gates, the Luttinger liquid is unsta- 
ble against the onset of a strong-coupling phase with very 
large CDW correlations. The temperature of crossover 
to the new phase depends crucially on the dielectric con- 
stant k of the environment, ranging from T c ~ 10~ 4 K 
(at k ~ 1) up to a value T c ~ 10 2 K (reached from 
K ~ 10). The inspection of the response functions mea- 
suring the superconducting correlations shows that these 
do not have an appreciable growth in the regime where 
the CDW correlations grow large. Thus, even under con- 
ditions of strong screening by external charges, the ob- 



servation of any superconducting feature seems to be ex- 
cluded in single nanotubes with the (5,0) geometry. 

The physical picture is different when we consider the 
case of a large array of nanotubes, of the kind patterned 
by the channels of the zeolite matrix in the experiments 
of Ref. 18. The coupling between the nanotubes in the 
array leads to a large reduction of the Coulomb potential, 
which is rendered finite at vanishing momentum-transfer. 
The Coulomb interaction happens to be now in the weak- 
coupling regime, making therefore easier the destabiliza- 
tion of the Luttinger liquid. We have found however that, 
before the superconducting correlations may grow large, 
a singularity is reached in one of the charge stiffnesses 
of the Luttinger liquid. This implies the divergence of 
the compressibility in one of the charge sectors, and it 
has therefore the same physical interpretation as the 
Wentzel-Bardeen singularity 23 . As the density-density 
correlator in each charge sector is proportional to the 
charge stiffness, the divergence of the latter leads to a 
phase where the system prefers to form electron aggre- 
gates, leaving some other regions with a defect of electron 
density. The possibility of a phase of this type has been 
already anticipated in Ref. [lT]- We have established that 
such a phase is realized under the particular experimen- 
tal conditions reported in Ref. 0, prevailing over any su- 
perconducting instability at the transition temperatures 
found in our analysis. 

The divergence of the charge stiffness opens the pos- 
sibility of having a genuine phase transition in the (5,0) 
nanotubes, since the above breakdown of the Luttinger 
liquid takes place without the formation of long-range 
order in a continuous order parameter. We have seen 
that, for a choice of the dielectric constant in accordance 
with the experimental conditions described in Ref. ITsl 
the transition temperature is of the same order as the 
temperature of the crossover to the pseudogap regime 
observed in that reference. It is worthwhile to stress that 
the divergence of the charge stiffness leads to the devel- 
opment of a pseudogap in the single-particle spectrum. 
The density of states near the Fermi level for electrons in 
subband i scales in the Luttinger liquid regime as 

n t (e)~e a > (49) 

with Ui as given in Appendix A. Quite remarkably, the I- 
V characteristics reported in Ref. 18, show the evolution 
of the curves for decreasing temperature into a power- 
law behavior with an increasingly large exponent. This 
is in agreement with the prediction obtained from 149|) 
when approaching the singularity in the charge stiffness. 
It becomes then plausible that the crossover to the pseu- 
dogap regime reported in Ref. may correspond to the 
breakdown of the Luttinger liquid found in the present 
framework. 

There are however other experimental features of the 
small-diameter nanotubes described in Ref. that can- 
not be explained in terms of the electronic instabilities we 
have found in the present paper. The strong diamagnetic 
signal measured at low temperatures has been interpreted 



12 



in that reference as a consequence of entering a super- 
conducting regime in the small-diameter nanotubes. It 
remains open the question of whether the armchair (3,3) 
nanotubes can be responsible for such a regime, or for 
some other type of electronic instability that may account 
for the strong diamagnetic signal, as some preliminary 
results seem to indicate^. Other questions that deserve 
attention refer to the properties of the small-diameter 
nanotubes when eventually assembled into other aggre- 
gates different to those considered in the present paper. 
It may be interesting to see the influence of intertube 
tunneling between small-diameter nanotubes packed to 
form massive ropes. It will be crucial then to discern 
whether the zigzag (5,0) and the armchair (3,3) nan- 
otubes lead to different phases, and whether they may 
support the appearance of larger superconducting cor- 
relations than those observed in massive ropes made of 
typical nanotubes. 
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The renormalized Fermi velocities of the new free-boson 
hamiltonian read: 



-1/2 
-1/2 



uq -> u /v- 



(50) 



Now we use this result to calculate the nonperturbative 
contribution A to the scaling equations of up and Ub- 
For instance, let us focus on up, whose operator has the 
form 



Up 



^ R0a ^R2^{ 0a ^Ll,'+h.C.). (51) 

The scaling dimension of H UF is readily evaluated with 
the bosonization technique, that allows us to express the 
fermion fields in terms of boson operators. For example 
let us bosonize one of the contributions to H UF : 



Appendix A: Diagonalization of Hfs and evaluation 
of the anomalous dimensions 

In this Appendix we diagonalize Hps by means of 
a generalized Bogoliubov transformation. Since in Eq. 
(|TT)l the {$_,II_} sector is already decoupled, it is suf- 
ficient to operate the transformation by mixing only 
the <E>+, <E>o fields and their corresponding conjugate mo- 
menta 11+ , IIo . Therefore we introduce the new operators 
$ + ,$o,n+,flo defined by: 
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where c (+) = cos<,so(+) and So(+) = sm( / J o(+) m 
order to ensure the standard commutation relations 
[$o {+) (x),fl(x') 0(+) ] =i8{x-x'). 

The angles ipo, ip+ and the parameters (j,, v, a, (3 depend 
on K , K + ,u , u+, f( 2 \f( 4 \ and must be determined by 
imposing that (i) [$(x) 0(+) , il(x') +{Q) } = 0, (ii) H FS 
written in terms of the new fields gets diagonal, (hi) the 
coefficients of (d x <& {+)) 2 and (n (+)) 2 are the same. The 
corresponding algebraic system has an analytical solution 
in closed form under the condition = /( 4 Ui: 



1 



-arctan 



3/2, 1/2 3/2, 1/2 



- I 1 



/u. 



& R0 ViuV] m *lz cx exp[--^(-0 flO + <f>m + 4>lo - 4>L2)} 

Now we can read the anomalous dimension of this oper- 
ator produced by the forward scattering interactions by 
expressing the boson fields <f> in terms of the operators 
which diagonalize Hps- The connection of the <j)'s with 
the operators entering in the hamiltonian in Eq. I|ll|) is 
given by: 
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faa = i(*+-*- + 0+-0-) 
0zi = ^(*+-*_-fl+ + e_) 

<t>R0 = ^($0 + #o) 

0LO = -($o-6 ), 
where the 9 fields are related to the IPs by: 

n±(o) = —dx8±(Q) ■ 

This allows us to write 

*L>**i*L>*£2 oc exp[-J=(0 o - -±=9+ - -j=B-)] 
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Finally the anomalous dimension A of the operator H UF 
is obtained by looking at the coefficients of the free fields 
in the exponent: 
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By means of the same reasoning, we can compute also 
the anomalous dimensions of the operators Oj^'Q defining 
the two-particle correlation functions, introduced in Ap- 
pendix B. They encode the nonperturbative corrections 
to the scaling Eqs. f53^l for the response functions x an d 
read: 



Appendix B: Scaling equations for the response 
functions 

In order to study the competition between the possible 
instabilities, one has to consider the scaling flow of the 
two-particle correlation functions %. They measure the 
strength of the quantum fluctuations for a given type of 
ordering induced by the pair-field , in such a way 
that 
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where L is the length of the nanotube, to m is the bosonic 
Matsubara frequency and P, Q = ±1 take into account 
the band-entanglement 32 . By neglecting the dependence 
on u> m and expressing the momenta according to Fig. (|14ll . 
the complete collection of the Fourier transforms of the 
response functions reads: 
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Finally we evaluate the anomalous dimensions gov- 
erning the density of states near the Fermi level for elec- 



P^ ma {p-k)a a /^ LW {p) 

= ^E [*U(P-*)^*^(P) 

k F R1 ) = 

^ E [*k«(p-*K ,/s *io/»(p) 

P,ck,/3 

+p^ R0a (p- k K' l3 ^(p)\ , 

Os$A k ~ °) = 7= E + fc)< /J *L2/3(p) 

P,Q!,/3 

*fl.2 Q (-p + fc)<^*Ll/3(p) 
+Q*flO«(p-fc)<' /3 *LO^(p)] , 

Osc>,v(k ~ k Frtl + k FL1 ) = 

■4 E Nat-p+^tivW] , 
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FIG. 14: Sketch of the Brillouin Zone of the (5,0) nanotube: 
the six Fermi points are indicated with the labels of corre- 
sponding chirality and subband. The longitudinal momentum 
transfers 2/cf and 2k2 are also drawn. 



Osc".^(k « k Flil + kp Ln ) = 



where for density wave (DW) operators /i = stands for 
CDW and ju — 1, 2, 3 for SDW; while for superconducting 
(SC) operators [i — stands for SS and fx = 1,2,3 for 
TS; cr^'' 3 are the Pauli matrices, with a^ 13 = 12x2- 

The response functions x defined in Eq. I|52|) do not 
obey scaling equations^, but the auxiliary functions de- 
fined as x(k, I) = TTVp(d/dl)x{k, I) do. Finally the scaling 
equations for the x functions read: 
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where the nonperturbative 
encoded in the anomalous 



AdW , AdW , AdW" , AnW" 

evaluated in Appendix A. 



Aa,b) 



, Asc, Asc 1 " 



(53) 

contributions 
dimensions 
are 



Appendix C: Intertube screening for forward 
scattering processes 

In this Appendix we study the Coulomb interaction 
at small momentum-transfer between electrons belong- 
ing to different nanotubes in a 3D array. The long-range 
intertube effects operate in two ways: i) it provides a 
screening of the bare intratube forward-scattering cou- 
plings gf> , , g { p , , g^ , g ( 4 ) , and ii) pro- 
duces new intertube couplings g^\g± (see Eqs. 1(15)1 . 
(|46|l ) with nontrivial scaling flow. Here we take into ac- 
count the screening effects by the nanotube environment 
by generalizing the approach devised in Ref. l33l which 
essentially consists on a RPA treatment of the dielectric 
constant. The RPA approximation is justified as long 
as the electrons are not allowed to tunnel between differ- 
ent nanotubes, and we assume that this condition applies 
to the system under consideration in the main body of 
the paper. It is important to stress that the screened 
interactions obtained by this approach do not contain 
any dependence on the energy cutoff. As a consequence, 
the RPA screening only imposes the initial value of the 
forward-scattering couplings in Section llVl but does not 
affect their scaling flow. 

Let 1 and 1' be the positions of two nanotubes mea- 
sured over a section of the 3D array, orthogonal to the 
tube axes. We introduce the long-range bare Coulomb 
interaction with momentum-transfer k along the longitu- 
dinal direction by partially Fourier-transforming the 3D 
Coulomb potential: 



Vi,v(k) = —K (k\\ - l'|), 

K 



(54) 



where K (x) is the modified Bessel function, which di- 
verges logarithmically as x — > and is exponentially 
suppressed for x > 1. It is worth to observe that V\y{k) 
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does not take into account the variation of the electronic 
Bloch wavefunctions around the waist of the nanotubes. 
This approximation is justified as long as the transverse 
momentum-transfer is small, i.e. for forward-scattering 
processes. When V\ t \'(k) is evaluated within the same 
nanotube, it is understood that |1 — l'| ~ R, and we 
recover the intratube interaction of Eq. with the 

logarithmic divergence at small momentum-transfer. 
The RPA-screened Vy[)(k) obeys the Dyson equation 

V$ (k) = Vw (k) + n £ V ir (k) V$\, (k), (55) 



where 1" runs over all the positions of the nanotubes in 
the 3D array and II is the polarizability of the ID electron 
gas; its zero-frequency Fourier transform reads 



n(fc) 



L 



E 



/(gq+fc) - f(£q) 
£q+k ~ £q 



(56) 



where / is the Fermi-Dirac distribution and L is the 
lcnght of the nanotubes. In order to solve Eq. I)55|l. 
we introduce the 2D Fourier transforms of V\ v (k) and 

V{;)(k): 



Vi,v(k) 



f d 2 p^(fc,p)e l P( 1 - 1 '), 

JBZ 



(57) 



and the same with V -> and 4> -> cj> {r \ where BZ 
indicates the first Brillouin zone of the nanotube lattice 
in the cross-section of the array, with lattice constant 



d ~ 1 nm. In terms of <p(k,p) and (/>( r )(fc, p), the Dyson 
equation reduces to 

(r) (fc,p) =#A,p)+n(fc)#fc >P )0M(*,p). (58) 

This permits to calculate the screened forward-scattering 
interaction with longitudinal momentum-transfer k be- 
tween electrons in tubes at 1 and 1': 



d 2 p 



<t>(k,p) 



BZ 



i-n(fc)<Kfc,p) 



ip.(i-io 



(59) 

The numerical evaluation of V^) (k) shows that such 
interaction is not negligible only for 1 = 1' and for nearest- 
neighbor tubes, i.e. for |1 — l'| = d. This calls for the 
introduction of intertube couplings with |1 — l'| = d. As 
discussed in Section ITVl two of them (g^ and g^p) have 
nontrivial scaling flow, and affect also the scaling of the 
intratube couplings. 

The logarithmic divergence of Vi t v (k) at k — is cured 
by the RPA screening; nevertheless the effective interac- 
tions are still enhanced at small k, but they saturate at a 
finite value. For the intratube interaction we have found 
that the saturation value is essentially independent on k 
and is V^\\k = 0)/vf ~ 0.65. In the intertube case, the 
RPA screening is much more efficent and the saturation 
value for |1-1'| = d is V$(k = Q)/v F ~ 0.002k, with an 
approximate linear dependence on k, at least up to k ~ 5. 
As discussed above, these asymptotic values are used as 
initial conditions for the forward-scattering couplings in 
the scaling equations in Section llVl 
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